Dissipation-Driven Breakdown of Universality in Two-Dimensional Superconductors 
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The influence of gapless dissipative degrees of freedom on the superconductor-insulator transition 
in two dimensions is investigated. We develop a series expansion for the free energy of a (2+1)- 
dimensional XY model coupled to a bosonic heat bath that can be approximately summed to all 
orders. The calculation explicitly conserves topological excitations. We derive the zero temperature 
phase diagram and the free energy critical exponent, and find a transition from universal to non- 
universal scaling behavior as the coupling to the dissipative environment is increased, implying the 
existence of a new universality class. 
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The subject of universality in the zero-temperature 
superconductor-insulator transition in two dimensions 
has been a debated issue for some time. Early exper- 
iments on thin superconducting films |Q seemed to in- 
dicate that the critical resistance right at the transi- 
tion is always in the vicinity of the resistance quantum 
Rq = h/(2e) 2 ~ 6.5fcil. This universality could be ex- 
plained by a scaling analysis of the bosonic Hubbard 
model |2|, which can be approximately mapped to the 
XY Hamiltonian, Eq. (|j). In later experiments very dif- 
ferent critical resistivities were measured in amorphous 
films U and Josephson junction arrays B, so that the 
universality hypothesis was cast into doubt. 

The universality class of the superconductor-insulator 
transition is affected by disorder and by coupling to dis- 
sipative degrees of freedom of the environment. In this 
paper we will concentrate on the latter, leaving the disor- 
dered model for future study. In a Josephson array with 
zero or weak dissipation the phase transition is driven 
by a competition between the Josephson coupling en- 
ergy and the charging energy and (at T = 0) belongs 
to the (2+l)d XY universality class. It is well known, 
on the other hand, that even in a single Josephson junc- 
tion the order parameter phase can localize if the elec- 
tronic degrees of freedom are taken into account, and if 
the associated density of states extends down to zero fre- 
quency 0. Clearly, under the same conditions, an array 
of such Josephson junctions would display phase order, 
and hence superconductivity, as well. Then the transition 
would have an essentially (0+l)-dimensional character, 
thereby motivating the existence of a new universality 
class. 

In this letter we will develop a framework in which 
these universality classes can be understood in a rigor- 
ous, quantitative way. Our formalism, which constitutes 
an expansion about the insulating state, explicitly pre- 
serves the 2n periodicity of the phase variables, and hence 
the topological excitations (vortices), which are known 
to drive the phase transition in the classical limit. We 
derive a series expansion for the free energy of the as- 
sociated (2+l)-dimensional statistical mechanics model, 




FIG. 1. The zero-temperature phase diagram in the 
non-intersecting path approximation. The X's in the inset 
show the free energy critical exponent as calculated here. The 
solid and dotted lines in the inset show the functional form 
given in Eq. hoi 



which can be approximately summed to all orders. At 
T = 0, this is of course just the ground state energy 
of the quantum system. Strikingly, we find a crossover 
from universal to non-universal scaling behavior as the 
coupling strength a to the environment is increased, as 
shown in the inset of Fig. (|l|): As long as a < 2/3, the 
free energy critical exponent is unchanged from its value 
for the non-dissipative model, while it varies continuously 
for a > 2/3, implying a new universality class. Fig. (||) 
shows a conjectured renormalization group (RG) flow di- 
agram which we will comment on at the end of the paper. 

The model under consideration has been investigated 
in the past, often by approximations that, in one way 
or the other, amount to some variant of a self-consistent 
field approximation, or to linearizing the Josephson cou- 
pling term, and hence destroying the 27T periodicity. 
However, a recent investigation [Q utilizing a Ginzburg- 
Landau- Wilson formulation showed results very similar 
to ours. 
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A number of mechanisms (see e.g. M for references) 
can give rise to an ungapped electronic density of states, 
and hence Ohmic damping, in thin-film superconductors. 
As a scenario that recently gained experimental support 
H we mention the formation of local pools of unpaired 
electrons due to spatial fluctuations of the order param- 
eter amplitude, caused by impurities. Other possibilities 
include d-wave high-temperature superconductors with 
nodal order parameter, coupling to electronic degrees of 
freedom in the substrate and in the vortex cores, or An- 
dreev scattering. A Josephson array with tunable cou- 
pling to the environment in the form of a 2d electron gas 
has actually been constructed by Clarke et al. ||] 

Two-dimensional Josephson junction arrays, as well as 
superconducting thin films can be modeled as a system 
of coupled quantum rotors with the Hamiltonian (in the 
absence of dissipation) 
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where the phase <f> of the superconducting order param- 
eter is a compact variable on the interval [0, 27r], and the 
Cooper pair number operator n is its canonical conjugate. 
Ec = (2e) 2 /C and Ej are the charging energy and the 
Josephson coupling energy, respectively. In the above 
Hamiltonian the ground state is limited to integer fill- 
ing factors, and more generally we would have to replace 
n x by (n x — S) 2 , where S is some non-integer number. 
This difference is critical in the non-dissipative model [^) , 
where a non-integer filling factor increases charge fluctu- 
ations and hence enhances phase ordering. However, we 
will see below that this dependence on the filling factor 
disappears as soon as we include Ohmic dissipation. 

In an imaginary time path integral formalism the ac- 
tion for the dissipative system reads 
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where the last term arises from an integration over 
the dissipative degrees of freedom flEo) and represents 
Cooper pair breaking processes due to the interaction 
with the environment. An alternative model would cou- 
ple the dissipative degrees of freedom to the phase dif- 
ference across the junction, reflecting the existence of a 
normal conducting channel between the islands, but will 
not be considered here. For Ohmic damping the kernel 
D(t) is given by its Fourier transform 
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for \uj\ less than a high-energy cutoff uj c , which can be ap- 
proximately set equal to some physical cutoff, such as the 




FIG. 2. The conjectured renormalization group flow di- 
agram. The isolated star marks the critical point of the 
non-dissipative XY transition, the other stars label a criti- 
cal line. The diamonds mark a line of stable fixpoints at 
E.,/Ec = 0. 



bandwidth of localized electron pools. The dimensionless 
parameter a controls the strength of the dissipative cou- 
pling. 

In general the compactness of the angular variables 
causes <j> and <fi + 27r to be identified, so that <p T is peri- 
odic only up to an integer multiple of 2ir. We therefore 
define <f> XT = <f> XT — 2-k(3 m x T l where xr satisfy peri- 
odic boundary conditions and the m x are winding num- 
bers. An infrared divergence in the dissipative term of 
(^), however, effectively suppresses all nonzero values of 
m x . This leads to another simplification as well: The 
non-integer filling factors mentioned above give rise to 
a term iS ji 3 drij) XT = 2iri5m x in the action. Hence, as 
long as winding numbers are suppressed, non-integer fill- 
ing factors are of no consequence. Intuitively speaking, 
since in our model the number of Cooper pairs is not 
conserved, any dependence on the filling factor is washed 
out. 

We now express the partition function as a power series 
in Ej: 
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Here [xtyi] indicates a directed bond between two 
neighboring lattice sites, and the sum runs over all such 



directed bonds. 



Vx 



i) can 



be interpreted as a "charge density," with Xi and y,; in- 
dicating the position of positive and negative charges, 
respectively. Each directed pair [xiyi\ then constitutes a 
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"dipole." Ej/2 plays the role of a fugacity, and So is the 
action (||) with Ej = 0. We may call the resulting model 
the Coulomb dipole model, in analogy to the Coulomb 
gas model for the resistively shunted Josephson junction 
H . After integrating out the phase variables the action 
reads 
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with 



v(t) 
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where we have subtracted out an infinite constant that 
enforces "charge neutrality" for each lattice site. For the 
rest of the paper we will assume Ec <C w c . Then the 
interaction potential is well approximated by 
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where 7 = 2/ a. This amounts to replacing the kinetic 
energy term E7}lo 2 by an effective cutoff uj* = Ec/2j. 
Each term in (^) can be represented by diagrams such 
as those in Eq. ([?]). Here each arrow corresponds to 
a "dipole" as defined above, or in physical terms to a 
Cooper pair transfer event. Double arrows, such as those 
in D 2 , indicate two dipoles of opposite orientation occu- 
pying the same bond. Similarly, the first diagram in D 4c 
has four dipoles occupying the same bond. As in stan- 
dard diagrammatic perturbation theory we can use the 
linked cluster theorem to write the perturbation series 
for the free energy as a sum over connected diagrams. 
The diagrams for the first few terms in the expansion 
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follow directly from Eqs. ^ to (||): 
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Here 
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and fij is a shorthand for f(ji —Tj). The subtractions in 
D 4 f, and D 4c are corrections for double counting and are 
distributed such as to render the diagrams finite in the 
zero-temperature limit, for a < 2. Beyond this value each 
diagram is divergent at T — 0, reflecting the phase local- 
ization transition in a single junction mentioned above. 

While for n > 4 the diagrams are generally very hard 
to evaluate, there exists a subset which can be computed 
exactly (within numerical limits) to all orders. The fol- 
lowing results are based on this partial summation, and 
we will comment on the validity of this approximation be- 
low. The diagrams we consider are the non-intersecting 
loops such as D 4a in Eq. (Q). We then have 
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where (at zero temperature) 
g(u) = / dr-, 



lu*t\ + l] 7 
-T (1 — 7, iuj) + c.c. 
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The numerical prefactor is equal to the number of 
translationally distinct self-avoiding polygons of length 
n on a square lattice, which is known |E| to behave 
asymptotically as A 2n ~ [i 2n , with the connective con- 
stant [i ~ 2.65. Hence the free energy is, in this approx- 
imation, 
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where G 2n = (uj*) 2n ^ 1 dw g 2n (cu) can easily be eval- 
uated numerically. This series constitutes an expansion 
about a disordered fixpoint of the RG flow diagram. The 
free energy is singular at the phase transition, and ac- 
cording to the fundamental properties of analytic func- 
tions this non-analyticity determines the radius of con- 
vergence. We can therefore derive the phase boundary 
from the asymptotic behavior of the coefficients in (pf): 
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Numerically we find that G n 1/n -»• (7- l)/2 for large 
n. This relation holds to an accuracy of more than five 
digits, and is presumably exact. Recalling uj* = Ec/2j 
and 7 = 2/ a we find that the phase boundary is at 
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The phase diagram is shown in Fig. (|l|). In order to 
assess the quality of the non-intersecting path approxi- 
mation it is important to realize that we did not omit 
any terms in the series expansion, but merely replaced 
the more complicated diagrams with products of simpler 
ones. For example, the connected piece of diagram 
is replaced by (-D2) 2 • To estimate the errors introduced 
by this approximation we calculated all intersecting dia- 
grams up to sixth order by Monte Carlo integration [jll) . 
We find that in fourth order these corrections change the 
value of Gn" by a factor between 1.1 (for a — > 0) and 
1.25 (at a = 2). For the sixth-order term these factors are 
1.12 (a — > 0) and 1.38 (a = 2). The second-order term 
is exact. Because of the relatively small ratios we expect 
our approximation to overestimate the critical value of 
Ej/Ec somewhat, but not to alter the qualitative fea- 
tures of the phase diagram. In particular, it cannot affect 
the universality classes. 

We now turn to the free energy critical exponent. Near 
the phase transition the singular part of the free energy 
must be of the form (1 — y) x , up to logarithmic correc- 
tions, where y = (Ej/Ej ) 2 and x is some unknown 
exponent. Expanding this expression in terms of y we 
get / ~ J2n c 2"Z/ n with 
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The leading asymptotic behavior of the coefficients is 
given by lnc2 n ~ — (x + 1) Inn. The corresponding terms 
in (Ist) behave as 
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where the first relation can be found in the mathematical 
literature jl2j (the exponent a s is conjectured to be 1/2), 
and the second relation has been determined numerically. 
We expect the numerical estimates to be correct at least 
to three digits precision. Combining these terms we find 
that the free energy critical exponent is 
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which is shown in the inset of Fig. ([!]). Due to the high 
quality of the fits shown we assume that this is an exact 
result with only minor corrections close to a = 2/3. The 
consistency between our results and those of Ref. M is 
emphasized by the scaling relation z = 2x — 3 between 
x and the dynamical exponent z calculated there. Since 
all correlation functions and thermodynamic properties 
of the system, as well as transport properties, can be de- 
rived from the free energy, this transition from universal 
to non-universal scaling behavior should be observable in 
those quantities as well. 



In the language of the renormalization group, the crit- 
ical point which controls the phase transition of the non- 
dissipative system has a finite basin of attraction. Since 
the flow equations are presumably analytic functions of 
the coupling parameters this implies the existence of an 
additional critical point on the phase boundary, at a = 
2/3. The continuously varying critical exponent further- 
more suggests a line of critical points beyond that value. 
The first-order RG equations [§ dEj/dl = (l-2/a)Ej 
and da/dl = (/ is a scale factor for the high-energy cut- 
off) determine the lower portion of the RG flow diagram, 
and in particular imply that Ej/Eq = constitutes a 
stable fixed line. Together these arguments lead to the 
RG flow diagram shown in Fig. (||). 
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